Gravity-analogy in one-dimensional ideal Fermi fluids 
and Burgers' equation 
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An hydrodynamic description of a one-dimensional flow of an ideal Fermi fluid is constructed 
from a semiclassical approximation. For an initially fully degenerate fluid, Euler and continuity 
hydrodynamic equations are dual to two uncoupled inviscid Burgers' equations. Yet the price 
for the initial simplicity of the description is paid by the complexity of non-linear instabilities 
towards possible turbulent evolutions. Nevertheless, it is shown that linear long-wavelength density 
perturbations on a stationary flow are generically stable. Consequently, linear sound obeys a wave 
equation with analogy to gravity. The results have applications for ultra-cold atomic gases. 



The behavior of a quantum held in the curved space- 
time of a black hole can be modeled by sound in a conver- 
gent fluid flow. For this, it is sufficient to have a local cor- 
respondence of equations: On one side the equations gov- 
erning the propagation of a quantum field near the event 
horizon; On the other side those governing the sound 
motion near the sonic horizon of a transonic flow. This 
correspondence has been first considered for a transsonic 
and irrotational inviscid-fluid flow by Unruh in his origi- 
nal work [l| . The remarkable consequence of such corre- 
spondence is Unruh prediction of thermal acoustic radi- 
ation emitted in transonic conditions. This prediction of 
acoustic Hawking radiation (sometimes referred as sonic 
Hawking radiation) suggested a way to test black hole 
evaporation in the laboratory. 

Acoustic Hawking radiation had attracted considerable 
interest and many aspects of it has been explored, such 
as the trans-Planckian problem just to cite one (see for 
instance and references therein). Moreover, gravita- 
tional analogues has been considered in different systems 
and in particular in specific quantum fluids, such as su- 
perfluids helium-4, helium-3, atomic Bosc-Einstein con- 
densates and one-dimensional (ID) Bose or Fermi de- 
generate fluid (for gravitational analogues see 0, 0| and 
references therein). 

A microscopic model of acoustic Hawking radiation 
that offers a way to visualize the source of the thermal ra- 
diation is desirable. A good quantum fluid candidate for 
such a model is the degenerate weakly-interacting atomic 
Bose gas where the large-scale fluctuations physics is 
mostly understood. This quantum fluid in a toroidal trap 
has indeed earlier suggested as a possible analogue of a 
sonic black hole [j| where to observe Hawking-like fluc- 
tuations effects. 

An even more fundamental model where the micro- 
scopic calculations can be done ab-initio from the many- 
particle Schroedinger equation is that based on a one- 
dimensional (ID) ideal Fermi-degenerate fluid @. This 
model is significantly more elementary than that based 
on the weakly-interacting atomic Bose gas in the corre- 
sponding one spatial dimension since both the complex 
physics of the last one and the inter-particle interaction 



are absent. Nevertheless, even if the ab initio calculations 
on the 1 D ideal Fermi fluid are simple to follow, the corre- 
sponding hydrodynamic description is more tricky. Such 
problematic is introduced in the present work. It is shown 
here that for an initially fully degenerate fluid, the Euler 
and continuity equations are dual to two uncoupled invis- 
cid Burgers' equations. This suggests the possibility of 
experimentally simulating a quantum version of Burgers' 
equation with ultra-cold Fermi atomic gases. Moreover, 
it is show here that linear sound propagating on a sta- 
tionary flow background is normally not evolving towards 
turbulent solutions in contrast to other solutions of the 
inviscid Burger equation. This guaranties the existence 
of solutions for the sound wave equation and establishes 
the analogy to gravity. 

Semiclassic equations of motion - - The hydrody- 
namic equations for a Fermi-degenerate system of non- 
interacting particles are constructed starting from the 
lowest possible order of a semiclassical approximation. 
That is of a classical dynamical system of non-interacting 
particles described by a phase-space distribution /, which 
evolves in a classical way according to the Liouville equa- 
tion, and where quantum mechanics and Fermi statistics 
enters only in the choice of the initial distribution. It 
is assumed that / varies in the spatial coordinate very 
slowly compared to the local mean inter-particle sepa- 
ration. We shall indicate by f(x,p,t) such semiclassical 
distribution of the fluid particles in phase-space as func- 
tion of time t, position x and particle momentum p. The 
particle density and the total local momentum (mass cur- 
rent) as function of spatial and time coordinates and are 
given by 



n(x,t) = J f(x,p,t)dp 
U(x,t) = / pf(x,p,t)dp 



(1) 

(2) 



respectively. The equation of motion for f(x,p, t) is given 
by the Liouville equation 



Lf + JLf^Ln 

dt dx dp 
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where H = H(x,p,t) = p 2 /2m + V ex t(x,t) is the Hamil- 
tonian of a single particle moving in an external po- 
tential V ext (x,t). The general solution of the Liouville 
equation © is formally obtained by the knowledge of 
the general solution of the trajectories x(x',p',t) and 
p(x',p',t) that are solution of the Hamilton equation of 
motion with initial condition {x',p'} at t — to. Then 
f(x',p',t) — f (x(x',p', — t),p(x',p', — t)) is a formal so- 
lution of the above Liouville equation. For instance, in 
the case of absence of external forces the solution is given 
by f(x,p,t) = fo(x — pt/m,p). For other possible ap- 
proaches to the formal solution of the Liouville equation 
see Reference Q. 

Inviscid Burgers ' equation — We limit here our con- 
sideration to the case of an initial fully degenerate ideal 
Fermi fluid. The quantum mechanical degeneracy of such 
a fluid is introduced via an initial condition /o with sharp 
left and right Fermi momentum Pq(x) and p R (x). This 
is given by 



- p\ 



(4) 



where h denotes the Planck constant. Since /o is then 
uniquely defined by the value of the left and right Fermi 
momentum Pq(x) and p R (x) we can look for the equa- 
tions of motions for these quantities. Indeed, a formal 
solution of the Liouville equation (|3|) with ((4]) as initial 
condition is given by the distribution 

f(x,p,t) = ±-0{ P -p L (x,t)}Qlp R (x,t) -p] (5) 

provides the left and right Fermi velocities Vl(x, t) = 
PL(x,t)/m and v R (x,t) = p R (x,t)/m satisfy both and 
independently the inviscid ID Burgers' equation Q 
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—v R (x,t) + v R (x > t)—v R (x,t) 
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v L (x,t) + v L (x,t)—v L (x,t) 

d 



d V ext {x,t) 

dx m 
d V ext (x,t) 
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dx m 

Note however that in order (|5j) to be a solutions of 
the Liouville equation <j3j> one has to take the velocity 
fields Vl(x, t) and v R (x,t) as a multi- valued function of 
x. This is indeed a known fact not restricted to the 
dissipation-less case Q • In the following we shall assume 
to be close to a stationary distribution with single- valued 
Fermi edges where infinitesimal perturbations (sound) 
usually remain single valued for an arbitrarily long time. 

Hydrodynamic equations — We shall assume in the fol- 
lowing that the left and right Fermi velocities vl (x, t) and 
v R (x,t) are single-valued. The hydrodynamic equations 
follows by taking the sum and difference of the left and 
right Burgers' equations ((6]) and ([7]). We obtain in this 
way the equation for the mean Fermi velocity 

v F (x,t) = [v R (x,t) -v L (x,t)]/2 (8) 



and the equation for the mean velocity 

v(x,t) = [v L (x,t) +v R (x,t)}/2 



(9) 



given by 
d 
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VF(x,t) 



dx 



d . , d 

dl vM = -Tx 



[v(x,t)vF(x,t)] 

v 2 {x,t) v 2 F (x,t) V ext {x,t) 



+ 



(10) 

(11) 



Since the fluid density ([T]) and the momentum density 
@ arc respectively given by 



2m 
~h 

2m 2 



v F (x,t) (12) 
v(x, £)vf(x, t) = mv(x, t)n(x, t) (13) 



n(x, t) 
n{x,t) 



we can recognized from (fTO)) and (fTT|) the continuity and 
Euler equation, respectively given by 

d d 
—n(x,t) = -— [v(x,t)n(x 7 t)) 
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dl 



v(x, t) 



d_ 

dx 



v 2 (x,t) 



w[n(x, t)] 



V ext (x,t) 



(14) 
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with an enthalpy w given by 
w(n) — 



8m 2 



(16) 



The local value of the speed of sound is as usual then 
fixed by the relationship c 2 = n(dw/dn) and it is given 
by the local value of the Fermi velocity 



c = vp{x, t) 



(17) 



Sound propagation and gravity wave equation — In the 
following we linearize the hydrodynamic equation and 
derive the wave equation for sound in a stationary fluid 
flow characterized by density no and velocity vq which are 
both assumed to be constant in time. We shall indicate 
by 



Sn(x, i) — n{x, t) — no(x) 
5v{x, t) — v{x, i) — vq(x) 



(18) 
(19) 



the density and velocity displacement from the stationary 
solution of flow, respectively. The linearized continuity 
and Euler equation read 

d d 

— 5n(x, t) = -— [5v(x, t) n (x) + v (x) Sn(x, t)} (20) 
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_Sv(x,t) = -- 



vq(x) 5v(x, t) + 



no (at) 



8n(x, t) 



(21) 



where Cq{x) = no(x) w'[no(x)} is the square of the local 
value of the speed of sound. After introducing the ve- 
locity potential $(x, f) = J v(x',t)dx' and the velocity 
potential displacement 



<P(x,t) 



5v(x' , t)dx 



(22) 



and some manipulation, we obtain finally the standard 
wave equation for 4>{x, t) by eliminating Sn(x,t) in the 
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linearized Euler equation (12 1 [) using the linearized conti- 
nuity equation (|20|) 



8 d . , 



n (x) 
c 2 (x) 



= iLb^iV^ (23) 

The above equation is the ID version of Unruh's wave 
equation [l[ , that is precisely the equation for a massless 
scalar field in a geometry with a metric given by 



ds 2 



n(x) 
c(x) 



| c 2 (x)-^(a;) dt 2 + 2v (x)dtdx - dx 2 } (24) 



Peculiarity of the acoustic-wave equation - The 
derivation in the previous subsection is a standard pro- 
cedure and it is valid not only for the specific case of 
a non-interacting Fermi fluid in one dimension. What is 
peculiar of the non- interacting Fermi liquid in one dimen- 
sion is the separation of left and right movers, which is 
explicit in the two uncoupled Burgers' equations and in 
their linearized counterparts. The consequence of this is 
that there is not scattering. This may sound a tautology 
since the particles of the fluid are assumed to not scat- 
ter and to be non-interacting, but it is worth to see the 
consequences at the level of the hydrodynamic equations. 
The wave equation (I23|) can be rewritten for the present 
case as 

1 
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v F (x) 
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dx 



The left and right movers stationary sound solution of 
equation (|25|) reads 



(f>L(x,t) 
<f>R{x,t) 



exp 



exp 



iu){ t + 



dy 



iuj \ t 



v F {y) - v {y) 
dy 



(26) 
(27) 



v F {y) + v (y). 

as one can easily check by explicitly replace them into 
equation (|2"5j) . (and) or into the linearized Burgers' equa- 
tions, which reads 

d d 
—5v L (x 7 1) + vl(x)—5v L (x, t) 

+ 6v L (x,t)^v° L (x)=0 (28) 

d d 
—Sv R (x,t) + v%(x)—5v R (x,t) 

+ Sv R (x,t)^-v° R (x)=0 (29) 

The reason of this absence of coupling is also related to 
the constant of the ratio of density and speed of sound 
that is given by 

n[x) 2m 
c(x) h 



(30) 



In fact, other fluids with gravity- analogy do not have 
the property that ((26)) and ([27]) are exact solutions of 
the linearized hydrodynamic equation but are actually 
an approximation of it consisting in neglecting an effec- 
tive potential proportional to spatial derivative of the 
so called conformal factor £(x)n(x) (£ = h/mc being 
the vacuum correlation length). This effective potential 
can scatter sound and it is often neglected in discussing 
acoustic Hawking radiation. The conformal factor for a 
ID ideal Fermi fluid is £(x)n(x) = l/n. 

Peculiarity of stationary flow — A different but mi- 
nor peculiarity of the fluid in consideration concerns the 
stationary flow equations. Let us indicated by / the 
conserved particle current no(x)vo(x). Then stationary 
equations are given by 



Vq(x) + v F (x) 



2(/i - Vest) 



III 



v (x)v F (x) = —I 
2m 



(31) 
(32) 



where /i is an integration constant. Note the symmetry 
in which the flow velocity and the Fermi velocity appear 
in the above equations. Expressing the above equations 
using the left and right Fermi velocities 



Vr(x) = v F (x) + v {x) 
v L (x) = v F (x) - v (x) 



(33) 
(34) 



which represent the local edges of the velocity distribu- 
tion, we obtain two uncouple equation 



mv 2 R (x) = (1r 



m v\{x) 



Vext (x) 
Vext{x) 



(35) 
(36) 



These equations are of course just the energy conserva- 
tion for the left and right Fermi edge, where the left and 
right chemical potential are given respectively by 



(37) 
(38) 



The metric written in terms of the left and right Fermi 
velocities is given by 



ds 2 



2m 
~h 



vr vl dt 2 + (vr + vl) dtdx — dx 2 



(39) 



Event horizon — Now we consider the particular case 
of a sonic event horizon. Let be vl the downstream mode 
that has an horizon in x = 0. This is realized by taking 
V ex t(x) = Vmax — \mn 2 x 2 as external potential and [il = 
Vmax- By inserting these in equation (|36|) we obtain the 
equation for the left Fermi velocity in presence of the 
horizon 



d , . d [ I 2 . 1 o 2 
M VLM+ dx-\2 VLM -2 KX 



= (40) 
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A stationary solution for the left moving edge which is 
differentiable is given by 

4 0) = kx (41) 

The equation for a displacement of the left Fermi point 
Vl = 5vl + is given by 

d d 

—Sv L (x, t)+g^ [nxSv L (x, t)} = (42) 
and has the solution given by 

5vL(x,t) —i — exp[— iuit + i(cj//t) \n.(x/a)] (43) 

KX 

Breaking time and quantum diffusion — The inviscid 
Burgers' equation is often used to exemplify the forma- 
tion of shocks and of multivalued solutions. The most 
simple example of this is that of an initial plane wave 
5v(x, 0) = Svosm(qx) in the absence of external forces. 
The solution with this single- valued initial boundary con- 
dition evolves into a multivalued solution after the break- 
ing time T c = 1/qSvo. To observe such non-linear effect 
in the ID ideal Fermi fluid T c should at least be shorter 
that the " lifetime" associated with the linewidth of the el- 
ementary excitation which is of order of r = 2m /Hq. So 
it is natural to ask how much such mode with wavevec- 
tor q should be populated in order that the non-linear 
effect starts to be comparable with the lifetime r. The 
answer is about one atom per half a wavelength. In fact, 

A /2 

indicating AiV = J Q Sn(x)dx = 2m5vo/hq results that 
T c = 2m/hq 2 AN = 2ttt/AN. Like viscosity in Burgers' 
equation tends to counteract the formation of shocks so 
acts the quantum mechanical uncertainty between posi- 
tion and momentum in the ID ideal Fermi fluid. This 
reminds us the viscosity bound conjecture of P. K. Kov- 
tun, D. T. Son, and A. O. Starinets |9|. 



Conclusion — It is shown here that ID ideal Fermi flu- 
ids in smooth external potentials admit a linearized hy- 
drodynamic description for sound with a wave equation 
with analogy to gravity. Thus ultra-cold atomic Fermi 
degenerate gases of negligible interparticle interaction in 
elongated traps with high radial confinement can be used 
as well as other quantum ID fluids to test in the labora- 
tory the acoustic analogue of Hawking radiation. More- 
over, a typical non-linear effect in hydrodynamics, i.e. 
the initial evolution towards the formation of shocks may 
be experimentally studied. 
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